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We present a model for inclusive electron and neutrino scattering off nuclei paying special attention 
to the influence of in-medium effects on the quasi-elastic scattering and pion-production mechanisms. 
Our results for electron scattering off are compared to experimental data at beam energies 
ranging from 0.7 — f .5 GeV. The good description of electron scattering serves as a benchmark for 
neutrino scattering. 



I. INTRODUCTION 

A wealth of information on nucleon structure func- 
tions, resonance widths and masses has been gathered 
through electron scattering processes at intermediate en- 
ergies (~ 0.5 - 2 GeV) (cf. , e.g., the reviews [11,11). At 
these energies, the lepton-nucleon reaction is dominated 
by quasi-elastic scattering and the excitation of baryon 
resonances and meson production. Furthermore, inclu- 
sive scattering experiments with nuclear targets have 
been performed by several groups, for a recent review 
see Ref. Q. In such experiments mesons and resonances 
are excited inside the nuclear medium. There, resonances 
and nucleons acquire an additional complex self energy 
due to rescattering effects and correlations which leads to 
modified spectral functions. By direct comparison to the 
elementary scattering case, one expects to deduce such 
in- medium modifications. 

Probing nuclei with weak probes, e.g. neutrinos, al- 
lows to study, compared to purely electromagnetic in- 
teractions, even more aspects of nuclear and hadronic 
physics, e.g., axial form factors and the strange quark 
content of the nucleon. However, the primary aim of 
investigating the vA process is related to the interpreta- 
tion of present-day neutrino oscillation experiments. A 
quantitative understanding of the influence of nuclear ef- 
fects on the cross section is crucial since most of the ex- 
periments use nuclei as targets Q. In this respect, the 
description of electron induced processes can serve as a 
benchmark for the neutrino induced reactions. 

There are three key issues in the theoretical un- 
derstanding of lepton-scattering off nuclei. First and 
foremost, one needs to model effectively the nuclear 
ground state and take into account the modification of 
the elementary lepton-nucleon vertex within the nuclear 
medium. Furthermore, the study of exclusive channels 
such as pion production or nucleon knockout demands 
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for a proper description of final state interactions of 
the produced particles with the nuclear medium. The 
latter issue can, e.g., be addressed within our Giessen 
BUU (GiBUU) framework [1, i, 0, B @ • 

There is a considerable amount of theoretical work 
aiming at a good description of the inclusive electron and 
neutrino cross section. Benhar et al. employ the im- 
pulse approximation with realistic spectral functions ob- 
tained from electron-induced proton knockout data and 
theoretical calculations based on nuclear many body the- 
ory (NMBT). With this model, they achieve impressive 
agreement in the quasi-elastic (QE) peak region; how- 
ever, they underestimate the data in the A region. In 
Refs. [nl, Ell they improved on this and a good descrip- 
tion of the data also in the single-pion production region 
could be reached. Also, Szczerbinska et al. use Ben- 
har's spectral functions for the QE contribution, but 
in the A region they apply the dynamical Sato-Lee model 
developed to describe photo- and electron-induced pion- 
production off the nucleon. Recently, this model has been 
extended to weak-interaction processes [H,!!!. A differ- 
ent approach to a combined investigation of neutrino and 
electron interactions makes use of the superscaling prop- 
erties of the electron scattering data (cf. Ref. [l6| and 
references therein) . There the authors extract the scaling 
function from inclusive electron-nucleus scattering data 
and use this to predict the neutrino-nucleus cross sec- 
tions. More work has been done in the QE region. Nieves 
et al. [l3| extended the nuclear inclusive electron scatter- 
ing model of Gil et al. [Tsll to electroweak probes which 
is in particular successful in describing the dip region in 
between the QE and A peak. Also Meucci and collabora- 
tors [H, HI apply a model developed for inclusive (e,e') 
reactions — a relativistic Green's function approach — 
to QE neutrino-nucleus scattering. 

In Ref. ^ the inclusive charged current neutrino cross 
section for both quasi-elastic and A region was calcu- 
lated, taking into account medium effects like the col- 
lisional broadening of the A. It is a mandatory check 
for this model that also total, inclusive cross sections 
for electron-induced reactions are well described. Thus, 
the present article aims at a combined study of inclu- 
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sive electron- and neutrino-scattering off nuclei, including 
both the QE peak and single-pion production. First, we 
present our description of the lepton-nucleus interaction 
in impulse approximation. Hereafter, we introduce our 
model for the nuclear ground state, for the in-medium 
modifications of the nucleon and the resonances. We 
discuss the influence of these in-medium effects on the 
lepton-nucleus cross section. We then present our results 
for electron- and neutrino- induced inclusive cross sections 
and compare these to available experimental data. First 
electron results have already been presented in plj . 



the appearance of effective masses for the nucleons. We 
include their effect both in the incoming current as well 
as in the final-state phase space. In the above equation, 
M is the mass of the incoming nucleon and ^(£",p') 
gives the spectral function for the outgoing nucleon {A ~ 
(5(p'^ — M'^) for free nucleons). 

The spin summed and averaged matrix element 
squared |A^p is proportional to the contraction of lep- 
tonic L^^ and hadronic tensor _ff 



(4) 



II. LEPTON-NUCLEUS INTERACTION 

In this section we present our model for the scat- 
tering of leptons with nucleons embedded in a nuclear 
medium. In particular, we discuss the electromagnetic 
(EM) {l~N l-X) and charged current (CC) (i^N 
l~ X) reaction on nuclei taking into accomit in-medium 
modifications. 

We treat the nucleus as a local Fermi gas of nucleons. 
The total reaction rate is given by an incoherent sum 
over all nucleons (impulse approximation) 
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where the cross sections on the rhs are medium-modified 
(cf. next section for details). Here we use the following 
notation: a lepton with four- momentum k = (fco, k) scat- 
ters off a nucleon with momentum p = {E, p), going into 
a lepton with momentum k' = (fcQ,k'). We further de- 
fine the transferred energy w = fcg — fcg, the transferred 
four-momentum = —q^ = — (fc — /c')^ and the solid 
angle n = Z{k,k'). 

At the energy region of interest (fco ^0.5 — 2 GeV), 
the cross section is dominated by quasi-elastic (QE) scat- 
tering {eN — > e'N' and vN — > l~N' , respectively) and 
single-pion production (eA'' — * e'nN' and vN — > 1~'kN'). 
Thus we assume 
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As we shall see in the next section, at higher lepton en- 
ergies the single-pion contribution is not sufficient to de- 
scribe the data reasonable well. Therefore, we shall ex- 
tent our model in this respect HI] . 

The cross section for QE-EM/CC scattering is given 

by 



dw dQ 



(3) 

where p' = (E' , p') is the four- momentum of the outgoing 
nucleon; mi is the mass of the incoming lepton. The 
presence of a momentum-dependent mean field leads to 



where the couphng is given by Cem — {4:na)^/Q'^ for 
EM and Ccc = Gl,cos^0c/2 for CC reactions. For the 
QE hadronic tensor we refer the reader to our earlier 
work in Ref. Q keeping in mind, that, in the case of EM 
scattering, the axial parts should be omitted. The vec- 
tor form factors are taken from the analysis of Bradford 
et al. [11] (BBBA-2005 form factors) and a dipole ansatz 
with Ma = 1 GeV is used for the axial ones. 

The single pion cross section dcr^M.cc/'^'^^^ domi- 
nated by the excitation of the A resonance and its sub- 
sequent decay: IN — > I' A I'ttN'. However, contri- 
butions to single-pion production from higher resonances 
and non-resonant terms are non-negligible. Henceforth, 
we introduce a background cross section which absorbs 
the latter contributions as well as their interferences with 
the A excitation term. Thus, the total single-pion cross 
section reads 
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To obtain da^M cc/^''^^^^ ^^'^ the same formalism 
as for QE scattering (see Eq. ^) only replacing H^'^ 
in Eq. (|4]) by the corresponding one for A excitation 
(cf. Eq. (18) in our earlier work [9[). Again, we omit the 
axial parts in the case of EM scattering. This method al- 
lows us to account for in-medium effects such as effective 
masses and collisional broadening as we shall discuss in 
more detail in the next section. Following Refs. [2^ [25j. 
we can relate the vector form factors to helicity ampli- 
tudes for which we take the results of the recent MAID 
analysis f2^ while the axial form factors follow a sim- 
ple dipole ansatz with an axial coupling obtained from 
PCAC. 

For the EM background contribution da^l^/dcudn to 
the single-pion cross section we subtract the dominant A 
contribution from the total single-pion cross section 



da 



A 



duo dO. duj dfl duo dO. 



(6) 



and obtain in this way the missing background. The 
total single-pion production cross section on the nu- 
cleon da]^/ duodVl is extracted from data. The elemen- 
tary vertex is expressed in terms of invariant ampli- 
tudes (27j . for which we use a parametrization of the 
MAID grou p [28l [29j (more details can be found in our 
earlier work |2l|). 
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FIG. 1: Elementary cross section for e~p e~ X for a beam 
energy of 730 MeV and a scattering angle of 37.1° in lab co- 
ordinates. The full curve denotes the total Ivr-cross section 
according to MAID. The dashed curve gives the contribution 
of the A resonance alone; it includes a proper real part accord- 
ing to dispersion relations. Data are taken from O'Connell et 
al. [30I ]: the statistical errors of the data points are negligible 
and, therefore, not shown. 



This procedure is visualized in Fig. [TJ The dashed 
line shows the A contribution to the cross section, 
i.e. da^^/dujdQ. Adding the background to the calcu- 
lation in the way described above (solid curve), a very 
good agreement with the experimental data of O'Connell 
et al. [30| is obtained. We conclude, that the A contri- 
bution alone is clearly insufhcient to describe the data 
accurately and that a background is needed. 

Such a treatment is not possible in the neutrino case, 
since there are not enough experimental data to fix the 
additionally necessary six axial amplitudes, hence we set 
dcr^%/dujdQ to zero. However, we plan to improve on 
that in the future (2^ . Since a phenomenological ansatz 
is not possible for CC reactions, elementary models are 
required to estimate the single-pion contribution from 
higher resonances and non- resonant terms [l^, [Sll, H^l ■ 



III. IN-MEDIUM MODIFICATIONS 

The target nucleus is treated within a local Thomas- 
Fermi approximation as a Fermi gas of nucleons bound 
by a mean-field. The nucleon mean-fiel d p otential is 
parametrized according to Welke et al. [33| as a sum 
of a Skyrme term depending only on density and a 
momentum-dependent contribution. Its parameters are 
fitted to nuclear matter properties in Ref. [s^l- In this 
case, M in Eq. ([3]) denotes the effective mass of the in- 
coming nucleon N, defined as 



M = Mn + Us^ip,r) 



(7) 



where Mn denotes its vacuum mass and C/s„ (p,r) the 
scalar potential. The spectral function A{E' , p') includes 
the effect of the momentum-dependent potential on the 
outgoing nucleon and also accounts for the in-medium 



collisional broadening of the outgoing final states. We 
neglect the spectral functions of the initial states because 
their widths are considerably smaller than those of the 
outgoing nucleons. 

Phenomenology tells us that the A(P33(1232)) reso- 
nance potential has a depth of about —30 MeV at po [ssj . 
in contrast to a momentum independent nucleon poten- 
tial, which is approximately —50 MeV deep. Therefore 
it is taken to be Va(p,i") = 2/3 ^(p,r). 

The density profile for g^O is chosen according to the 
parametrizations collected in Ref. (36| . The proton den- 
sities are based on the compilation of data from elec- 
tron scattering [37| : the neutron densities are provided 
by Hartree-Fock calculations. 

The spectral function of a particle with four- 



momentum p 
by 

A{E,p) = - 



{E, p) and mass M 
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with the self energy S(i?, p) and the vacuum pole-mass 
Mq. The imaginary part of the self energy is determined 
by the full width, Ftot, in the medium 



ImI](P,p) - -MFtot iE,p) 



(9) 



To deduce this width, we have to consider on one side 
the modification of the free A decay- width rY-oo, which 
is parametrized according to Manley et al. [38|. Due to 
Pauli blocking of the final state particles in the medium, 
the free decay width is lowered and has to be replaced 
by the Pauli blocked decay width: Ffrec ^pb- On the 
other side, both the nucleons and the A resonances un- 
dergo collisions with the nucleons in the Fermi sea. This 
leads to a collisional broadening of the particle width. 
To estimate the collisional broadening, we employ the 
low-density approximation 



F, 



□ii(£',p) = J n{p) (t{E, p,p') p Urci PpB 



dV, (10) 



with the momentum distribution n(j}) of the nucleons in 
the Fermi sea. Furthermore, (t{E, p, p') denotes the total 
cross section for the scattering of the outgoing nucleon 
(or A) with a nucleon of momentum p' in the vacuum. 
The variable Vrei denotes the relative velocity of the parti- 
cle and the nucleon; and PpB is the Pauli blocking factor 
for the final state particles. The total cross sections are 
chosen according to the GiBUU collision term d- Alto- 
gether, the full width is given by 



Ft, 



PB 



F, 



(11) 



In the same line as Ref. [33] , we demand that the real part 
of the self-energy at the pole energy is given by the mean- 
fields. Off the pole, the real part of the self energy is 
extracted by a once-subtracted dispersion relation. This 
procedure guarantees the normalization of the spectral 
functions. 
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The amplitudes for resonance production and quasi- 
elastic scattering are evaluated with full in-medium kine- 
matics. Furthermore, also the flux and phase-space fac- 
tors are evaluated with in-medium four-vectors and the 
spectral functions of the outgoing nucleons. As an ap- 
proximation, we use in the medium the same form-factor 
parametrizations as in vacuum. Pauli blocking is taken 
into account by multiplying each cross section on the rhs 
of Eq. (m with the Pauli-blocking factor. 

For the single-7r background we do not apply any in- 
medium modifications besides Pauli blocking and Fermi 
motion. When evaluating Eq. ([6]), we assume vacuum 
kinematics for the in- and outgoing nuclcon and pion. 



IV. RESULTS 



In Fig. [21 we present our results for the inclusive reac- 
tion l^O (e~, e~) X. We compare them to data by Anghi- 
nolfi et al. [401 at beam energies ranging from 700 to 1500 
MeV and a fixed electron scattering angle of 6k' = 32°. 
The long-dashed curves in all three panels show our result 
with Fermi motion and Pauli blocking, but without any 
mean field. Especially at 700 MeV, the QE peak is over- 
estimated and the so-called "dip- region" in-between QE 
peak and A peak is significantly underestimated. When 
the momentum dependent mean field is included (short- 
dashed line), then the faster (on average) final state nu- 
cleons experience a shallower potential than the initial 
state nucleons. Therefore, more energy must be trans- 
ferred by the photon such that the energy is conserved 
in the QE reaction. Hence, the QE peak is broadened 
towards a higher energy transfer w. By a similar effect, 
the single- TT contribution is also shifted in energy, because 
the outgoing A experiences a less attractive potential. Fi- 
nally, the solid curve represents a full calculation which 
includes both mean field and in-medium spectral func- 
tions incorporating proper real parts for the nucleon and 
A self energies. The QE peak is broadened even more 
and, as a consequence, the peak height decreases. The 
modification of the A-width in conjunction with a proper 
normalization has, however, only a minor impact on the 
total cross sections. 

We conclude, that the overall agreement to the data 
is improved by a calculation which in addition to a local 
Fermi gas momentum distribution also includes a mean 
field and in-medium spectral functions for the nucleon. 
Especially at low energies (cf. the upper panel represent- 
ing 700 MeV beam energy), a proper treatment of the 
nucleon spectral function is important. The increase of 
the energy loss due to a momentum-dependent nucleon 
potential reshapes the QE peak considerably. An addi- 
tional modification of the nucleon width leads to further 
broadening and decrease of the QE peak height. In the 
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FIG. 2: In the upper figure, the inclusive electron cross section 
da-E.m/ dujdO, on ^''O is presented as a function of the energy 
transfer = fco — fep at a fixed electron energy and scattering 
angle of 9y = 32°. The long-dashed line denotes our result, 
where we include Fermi motion and Pauli blocking only. The 
short-dashed line denotes the result, where we in addition take 
into account the binding in a density and momentum depen- 
dent mean-field potential. The solid line, finally, includes in 
addition the in-medium spectral function (SF). The A con- 
tribution to the full calculation (solid line) is denoted by the 
dotted hne. The data are taken from Ref. [4Cll. 1411. 



single- TT and A production region^, we achieve a good de- 
scription for all energies. The in-medium modifications 
improve the overall correspondence with the data. In the 



^ In the medium, the A has also pion-less decay modes and con- 
tributes, therefore, not only to single- vr production. 
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here, that the high w-region is underestimated since 
higher-mass resonances and background contributions 
are still missing. For the chosen observable, no exper- 
imental data are available and a direct comparison in 
order to check the quality of our calculation is therefore 
not possible. We may, therefore, use our electron scat- 
tering results as benchmark, where the same model with 
the same assumptions gives a satisfying agreement to the 
data, and thus serves as a quality check also for the neu- 
trino induced reactions. 



FIG. 3: Same as Fig. [5] but for i/^ CC scattering on ^"^O, 
= 1 GeV and 6 = 30°. 



dip region, which is conventionally attributed to 2N exci- 
tations, the description is considerably improved due to 
the previously discussed broadening of the QE peak. At 
higher beam energies (1.08 GeV, cf. middle panel; 1.5 
GeV, cf. lower panel), the data are underestimated at 
high uj due to the fact that 27r-production channels have 
not yet been included. 

An impulse approximation calculation by Benhar 
et al. [lo| that uses NMBT spectral functions yields in 
the QE-region a better result for 700 MeV beam energy. 
However, already at a slightly higher beam energy of 1080 
MeV our model and the NMBT one yield equally good 
results for the QE peak. We thus conclude, that our sim- 
ple ansatz for the in- medium width (cf. Eq. PH)) ) and 
the inclusion of a proper potential incorporate the main 
features of the nucleon spectral function in the medium. 
The NMBT pion-contribution of has lately been im- 
proved in [ 12 1 using similar methods as in our calculation 
of Ref. IJ"^ 

In Fig. IHl we present results for the CC reaction 
1^0 (i^e, e~) X for a neutrino beam energy of 1 GeV and 
a lepton scattering angle of 30°. Here we have further re- 
fined the results for the inclusive cross section shown in 
Ref. by including in-medium spectral functions. The 
treatment of the in-medium effects is as that used for 
the electro-production calculations described above. In 
analogy to the electron case, we show results without 
(dashed) and with (short dashed) mean field and the 
full calculation which includes both mean field and in- 
medium spectral functions (solid curve). One observes 
similar features as in the case of electron scattering: 
A broadening and a shift of the QE and singlc-7r peak 
caused by the momentum-dependent potential and the 
in- medium width of the nucleon and the A. 

Comparing our QE result (left peak) with the QE cross 
section obtained by Benhar et al. (their Fig. 13) — 
they apply the same model as in the electron case — 
we find that our results are slightly above their calcula- 
tion, but otherwise reproduces the main features of their 
result. 

As in the case of electron scattering, we also expect 



V. SUMMARY AND OUTLOOK 

In this work, we have presented a model for inclusive 
electron and neutrino scattering off nuclei. In particu- 
lar, we have studied the influence of in-medium modifi- 
cations on the QE scattering and pion-production mech- 
anisms at intermediate lepton energies. QE scattering 
is described with a relativistic formalism that incorpo- 
rates recent form factor parametrizations. We take into 
account the full in-medium kinematics, mean-field poten- 
tials and the nucleon in-medium spectral function. For 
pion production in e~-induced reactions, we have used a 
treatment based upon the MAID analysis, which allows 
to incorporate background terms besides the dominant 
A resonance. Here we also consider in-medium kinemat- 
ics, potentials and the in-medium modification of the A 
spectral function. For neutrino-induced pion production, 
however, we assume pure A dominance and neglect back- 
ground terms. The vector N — transition form factors 
are based on the latest MAID analysis. 

For electron induced processes at beam energies rang- 
ing from 0.7 — 1.5 GeV, we achieve good agreement both 
in the QE region as well as in the pion production region. 
The inclusion of mean field potentials and in-medium 
spectral functions improves considerably the correspon- 
dence to the experimental data obtained by Anghinolfi 
et al. [13] ■ A comparison to the results in the QE re- 
gion obtained within a NMBT calculation by Benhar 
et al. [10] shows that our implementation of in-medium- 
modifications provides an effective and efficient treat- 
ment of nucleon properties in the medium. When the 
energy transfer to the nucleus becomes high enough, 27r 
production sets in. In that region, we underestimate the 
data due to a lack of those mechanisms in our model. 

Considering the electron results as a benchmark, we 
have presented a calculation at 1 GeV neutrino energy 
relevant for current and future neutrino oscillation ex- 
periments. 

Our GiBUU transport model [5j also allows to study 
exclusive reactions, like pion production and nucleon 
knockout, taking into account re-scattering effects lead- 
ing to a change in the final state particle multiplicities 
and distributions. Using this model, we have already 
performed detailed calculations for exclusive, neutrino- 
induced, reactions like pion production and nucleon 
knockout [1, Q. These calculation use almost exactly 
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the same treatment for the initial interaction process as 
described here. In the future, we plan to improve on the 
results of Lehr et al. [13] where exclusive channels such 
as 77, TT and tttt have been evaluated for electron-nucleus 
scattering within a precursor version of GiBUU. 
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